Kondo Scaling in the Optical Response of YbIni_j;Aga;Cu4 
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Theoretical work on Kondo systems predicts universality in the scaling of observable quantities 
with the Kondo temperature, Tk- Here we report infrared-frequency optical response measure- 
ments of the correlated system YbIni_2:Agi,Cu4. We observe that x-dependent variations in the 
frequency and strength of a low energy excitation are related to the x-dependent Kondo temper- 
ature. Comparison of the inferred trends with existing theory and a model calculation provides 
a framework in which to view these experimental results as scaling phenomena arising from local- 
moment /conduct ion electron hybridization. 



The study of the Kondo problem has contributed much 
in the way of theoretical technology (e.g. renormalization 
group) as well as insight into the possible manifestations 
and phenomenologies of interacting systems with many 
degrees of freedom. One unifying effect born from the 
years of investigation is the Kondo resonance-a many- 
body collective excitation responsible for much of the in- 
teresting low energy behavior displayed by rare-earth and 
transition metal systems. 

YbAgCu4 provides an example of a moderately heavy 
fermion system. At low temperatures, the j = \ Yb 
moment is screened by conduction electrons fflving an 
enhanced Pauli paramagnetic susceptibilityP, En- 
hancement is also found in the Sommerfeld coefficient 
7 ~ 250 TO J/mo^-ftT^, indicating a large effective carrier 
mass0, or equivalently, a large density of states near 
the chemical potential. At high temperature, the Yb mo- 
ment acts independently of the carriers, which show un- 
enhanced mass behavior. The crossover between the high 
and low temperature regimes occurs around a character- 
istic scale Tk, which is intimately related to the degree 
to which the low temperature properties are enhanced. 

Substituting In for Ag gives very different behavior. 
YbInCu4 displays a first-order isostructural electronic 
phase transition at T^, ~ 42 X from a free-moment 
semimetal to a low-temperature metallic mixed-valent 
phased 0. This phase transition begins a phase bound- 
ary in the x — T plane of YbIni_3;Aga;Cu4 which in- 
creases in temperature as x is increased and ends in a 
critical point around x ~ 0.3. This phase boundary is of 
considerable interest SB because it involves 
mainly electronic degrees of freedom and because there 
is as yet no real understanding of or agreement regard- 
ing its fundamental origin. When approaching an under- 
standing, one must address the periodic Anderson model 
(PAM), which provides a powerful basis for interpreting 



the phenomenology of heavy fermion and mixed-valent 
systems, where crossover behavior is the norm. While 
much of the phenomenology of YbIni_j;Aga;Cu4 is con- 
sistent with the PAM, the occurrence of an unexpected 
phase transition at low x could indicate the need for an 
additional interaction term in its minimal model Hamil- 
tonian. It is natural in this context to consider whether 
adding conduction electron interactions (c./. Giamarchi 
et al. 8\ and Freericks et a/.jli]) which are not generally 
included in the PAM, could provide competing influences 
necessary to understand the origin of this phase transi- 
tion. Further exploration in this area may help estab- 
lish a connection between the domain of Mott-Hubbard 
physics, where strong conduction electron interactions 
lead to a phase transition, and the moment-compensation 
physics of the Anderson and Kondo models. 

In low-temperature YbInCu4, where the Kondo scale is 
large, there appears a distinct feature in the mid-infrared 
conductivity that is not present in the high temperature 
phase ■ It has been suggested that similar mid-infrared 
features found in other heavy fermion materials |^ 0, 
lioj are a manifestation of local moment (Kondo lat- 
tice/periodic Anderson model) phenomena [ill 
[l3 | . Here, by studying the Ag concentration dependence 
of the dynamical conductivity in YbIni_a;Ag2;Cu4, we ef- 
fectively vary Tk within the low temperature phase and 
observe systematic changes in characteristics of the mid- 
infrared excitation. We then explore the scaling behavior 
of this excitation in the context of a simple local-moment 
hybridization picture, and establish both the identifica- 
tion and phenomenology of the optical signature of the 
Kondo resonance. 

We have measured the reflectivity of single crystals of 
a sequence of samples of YbIni_a;Aga;Cu4 with x — 0, 
0.3, 0.5, 0.75 and 1.0. Our measurements cover the fre- 
quency range from 5 meV to 6.2 eV with detailed temper- 
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FIG. 1: The real part of the low-temperature optical conduc- 
tivity of YbIni_a;Agi:Cu4 (T = 20 K) is shown for five values 
of X. The open triangles mark peak frequencies. The dotted 
curves refer to a PAM-based joint density-of-states calculation 
in the text and the dark vertical bars indicate a threshold fre- 
quency from that calculation. 



ature dependent data taken between bmeV and 2.8 eV. 
A Kramers-Kronig transform is applied to the measured 
reflectivity to determine the real and imaginary parts 
of the dynamical conductivity. For the purposes of the 
transform, Hagen-Rubens terminations are used below 
bmeV. At high frequency (above &eV) each reflectivity 
spectrum is extrapolated to a common value of 0.08 at 
bj — \h eV and then continued as a constant to 25 eV. Be- 
tween 25 eV and 50 eV an form is used for R{uj) and 
above 50 eV the free-electron form w"'* is assumed. In or- 
der to assure that our essential results below 1 eV do not 
depend on the high-frequency extrapolations, we have ex- 
perimented with a number of termination protocols in- 
cluding other common values and coalescence frequencies 
as well as constant extrapolations above 6.2 eV. These 
show convincingly that our results regarding trends in the 
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FIG. 2: a) Characteristic frequencies, lo/s (A) and ujth (a) are 
shown as a function of x (inset) and the Kondo temperature, 
Tk. b) Spectral weight is shown versus x (inset) and Tk- Tk 
is obtained from low-temperature magnetic susceptibility Q 
measurements. Spectral weight is calculated using a bare 
band mass m = 4me. Dotted lines show fits to the data based 
on model equations in the text. 



X and Tk dependence of cri(cij) below leV are not sig- 
nificantly influenced by any of the extrapolations above 
6.2 eV. 

In Figure^ the low temperature conductivity is shown 
for the five x values studied. A prominent feature 
centered around ^eV for x—0 (YbInCu4) changes only 
slightly as x is increased to 0.3, however further doping 
to a::=0.5 and 0.75 causes both a red shift and strength- 
ening of this feature. The triangles mark the frequency 
of the peak of cti (uj) and the vertical bars mark the lo- 
cation of a threshold frequency obtained from a model 
calculation described below. 

In Figure 2a these two characteristic frequencies, cor- 
responding respectively to the peak and threshold of this 
excitation, arc plotted both as a function of x and as a 
function of Tk . Values of Tk as a function of x are ob- 
tained from measurements of magnetic susceptibility 
The complex dependence on x becomes monotonic when 
plotted vs Tk. The dotted curves in Figure 2a refer to a 
theoretical scaling described below. 

In Figure 2b an integral of ai{Lu) in the vicinity of 
this excitation is shown as a function of x and Tk. This 
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spectral weight, or strength, is defined by the integral. 



n{Lu) 



2m 



<Ji {Lo')duj' . 



(1) 



where the limits of integration, 6meV and 0.5 eV, en- 
compass the J eV excitation and exclude a very small 
Drude peak at very low frequency. Neither integration 
limit is critical; in fact a lower limit of and upper limits 
between anywhere between QAeV and leV produce the 
same x dependence which makes one confident that the 
n versus x dependence shown here is an essential charac- 
teristic of the data, and independent of any of the details 
of the choices we have made in the analysis. 

Our discussion of the data centers on the PAM disper- 
sion relations illustrated in the inset of Fig. 3a which ad- 
dress the essential physics of hybridization of conduction 
electrons with local moments and the appearance of the 
Kondo resonance at Ep- At energies far from the chem- 
ical potential, the upper and lower bands, e+ and e~ , 
follow closely the unrenormalized free carrier dispersion. 
Near Ep, the Fermi surface opens (and the bands fiatten) 
to accommodate the /-electron weight projected up to 
the Fermi level as a result of hybridization. This reorgani- 
zation of the bands in the vicinity of the Fermi level is due 
to many-body interactions, the strength of which is char- 
acterized by the parameter V. The resultant narrow peak 
in the density of states, called the Kondo, or Abrikosov- 
Suhl, resonance is central to the understanding of heavy 
fcrmion and mixed- valent Dhenomenologv [l5l lid IT^ . 

The hybridization-induced splitting creates the possi- 
bility of vertical transitions from filled states below Ep 
to unoccupied levels above Ep, as illustrated by the ar- 
rows of Figure 3a. The threshold for these transitions 
occurs at a frequency oj — 2V, where V is the hybridiza- 
tion strength renormalized by the on-site /-electron re- 
pulsion. This energy scales with the Kondo temperature 
as[lll[ll[ll[ll[l3 



(2) 



where B is related to the conduction electron 
bandwidthfl^. At threshold, the nesting of the coher- 
ent (lower) and incoherent (upper) bands leads to a very 
high joint density of states for vertical transitions and 
hence a strong peak in the conductivity, which we show 
below can be related with the width of the Kondo reso- 
nance, Tk- 

To proceed with the analysis, we use the Kubo- 
Greenwood formulail9j : 
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where Ip^ £'| denotes the dipole matrix element connect- 
ing electronic bands £ and and JDOSt/i{Lo) is the 
corresponding joint density of states. Applying this for- 
mula to hybridizing quasiparticles as though they were 




FIG. 3: The inset shows the PAM dispersion relations 
with arrows indicating transitions made allowed by the ex- 
istence of the Kondo resonance, a) shows crpam(i^) calcu- 
lated from those transitions, b) shows a model conductiv- 
ity calculated using V — 41meV, e/ — 2.5 meV, Ep = leV, 
IPo- J'^kp = 4.673l"^ and A = O.UeV. 



electrons allows an exploration of the phenomena of the 
mid-infrared conductivity in the context of the PAM. In 
that case there are two relevant bands, e"*" and and 
our model approach leads to: 



dS 



lP+, 



(4) 



where |p_|_ _| is the matrix element for the one remaining 
term in the JDOSi^g' sum of Eq. ©. 

For each value of w > 2V, there are two contributions 
to the integral: one originating from levels inside the 
unrenormalized Fermi surface (A:<); the other from the 
states occupied as a result of renormalization, i.e., outside 
the unrenormalized Fermi surface (fc>). An example of 
two such transitions with the same lu are indicated by 
vertical arrows in the inset of Figure |3| 

Simplifying equation ^ by making an isotropic ap- 
proximation and the approximation of a constant |p^ | 

(discussed further below), one obtains 



CTpam(t<-') 



E 



47rfc2 
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k=k' 



Model spectra can then be calculated using the explicit 
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PAM dispersion relation [lH.ll^ . 

i ^ ^ ^ ^ , (6) 

where Ep is the Fermi level, V is the renormalized hy- 
bridization strength, introduced above, and e/ is the /- 
level position renormalized by on-site /-electron repul- 
sion which defines the scale of the low-energy physics and 
is commonly identified with the impurity Kondo temper- 
ature, Tk- 

An example of such a calculation for V = Al meV and 
Cf = 2.hmeV, is shown in Figure 3a. The onset for op- 
tical transitions appears as a cusp at 2V'|23|- While 
the strength of this cusp is an important feature asso- 
ciated with nesting, the extreme sharpness is an arti- 
fact of assuming infinite lifetimes for the quasiparticle 
states. We rectify this by convoluting this model result 
with a Lorentzian of constant width to post hoc emulate 
the finite-lifetime effects. The magnitude of the broaden- 
ing parameter, A, gives an estimate of the quasiparticle 
lifetime. The values used here of around 0.12 eT^, corre- 
spond to a lifetime of order t = /i/A ~ 6ps, consistent 
with pump-probe lifetime measurements of YbAgCu4 . 

In Fig. 3b we present this calculated Upamd^) to- 
gether with our measured ai{uj) for a;=0.75. Based on 
the quality of this fit, the PAM parameters V = 41 meV 
and Cf = 2.5 meV can be associated with this composi- 
tion. With V and e/ as adjustable parameters, we have 
similarly modeled the mid-infrared conductivity (Ti(w) 
for each x value as shown by the dashed curves in Fig- 
ure ^ The vertical bars mark the threshold frequencies 
ujth — 2V. These threshold frequencies are plotted versus 
X and Tk in Figure 2a. The dotted curve through those 
data shows the consistency between our experimental re- 
sult and the theoretically based scaling of Eq. |21 

With regard to the total strength^ one can obtain an 
analytic result in our approach by setting the two sphere 
areas in © equal to Airkp. This approximation avoids 
the effects of bare band structure details while including 
the influence of the strong cusp at 2V. Solving for the 
wavevectors fc> and /c< using the condition e+ — e~ = cj, 
one can obtain a model strength, which is essentially the 
shaded area of Fig. 3a, of: 



npam = 5- / Crpam{oj)dLJ (7) 

T^e J2V 

where c is defined by Tk — c e/ . This logarithmic scaling 
behavior is compared to the data in Figure 2b. 

It is important to realize that the simple arguments 
leading up to equation (|5J) are not strictly rigorous and 



in fact careful analysis reveals that a k-independent hy- 
bridization leads to a |p_|_ _| which depends on k, and 
therefore w, in a manner that will reduce the conductiv- 
ity at high frequency O^i^ expects that in a real sys- 
tem, transitions between Cu-In-Ag d orbital bands and 
Yb /-band states are likely to be allowed due to their rel- 
ative spatial displacement. Incorporation of this possibil- 
ity into theory can be done by allowing a dispersive Vj^., 
which would tend to restore some of the high frequency 
conductivity that is underestimated in our simplified ap- 
proach. A more sophisticated approach to n{TK) would 
be valuable in the ongoing effort toward a quantitative 
understanding of strongly correlated electrodynamics. 

Equations (9) and (2) do, however, qualitatively ex- 
plain the observed dependence of the frequency and 
strength on Tk- In addition these equations show a in- 
teresting relationship between the scaling of n(TK) and 
^th{TK)- Using the impurity model result for c 21], and 
fitting the measured n{TK) to Eq. @, one gets a value 
for the band parameter B = 0.36 eF. This value is de- 
termined solely by the extrapolation of the data to the 
Tk axis, and so does not depend on kp or the matrix 
elements present in the coefficient of the logarithm. One 
can independently obtain a value for B from fitting the 
Tk dependence of the excitation frequency (Figure 2a), 
LOth = 2V = 2v^IVS, which leads to B = 0.30 eV. The 
observation that these independently obtained values of 
B are comparable shows a non-trivial relationship be- 
tween the scaling behavior of strength and characteristic 
frequency in our optical data. Overall our results in- 
dicate that Tk scaling is present in the low temperature 
finite-frequency dynamics of YbIni_3;Aga;Cu4 and can be 
addressed in the context of local-moment models. Fur- 
ther work may be directed toward greater understanding 
of this low-T scaling behavior as well as unexplained tem- 
perature dependent behavior of Yblui-ajAg^: CU4. 
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